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20.  Abstract  (continued) 

station  site,  the  superposition  equations  defining  each  case  are 
solved  in  closed  form  to  yield  the  unknown  complex  source 
amplitudes  and  source  azimuths. 

Actual  filter  procedure  involves  dividing  the  component  signals 
into  time  intervals,  taking  the  Fourier  Transform  for  each  time 
interval,  and  operating  on  each  of  the  frequency  components  in  a 
time  interval  with  the  ellipticity  filter.  The  source  waveforms 
are  recreated  by  an  inverse  transform  of  those  complex  ampli- 
tudes associated  with  the  same  azimuth. 

Preliminary  evaluation  of  this  technique  has  been  performed 
and  in  the  noiseless  case,  has  perfectly  separated  two  mixed 
Rayleigh  waves. 
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1.0  Introduction  and  Summary 


One  of  the  important  aspects  of  any  underground  nuclear  test 
ban  treaty  is  the  ability  to  verify  that  the  signators  are  adhering  to  the 
agreed  provisions.  To  this  end,  reliable  methods  for  detecting  under- 
ground nuclear  events  must  be  employed.  A method  which  has  proven 
to  be  quite  useful  in  this  application  has  been  seismic  verification. 

Seismic  verification  of  a test  ban  on  underground  nuclear 
explosions  requires  the  classification  of  seismic  events  at  teleseismic 
(large)  distances.  Classification  is  determined  by  examining  properties 
of  the  seismic  waves  as  recorded  at  seismometer  stations  located 
around  the  world.  Two  wave  types  are  important  in  the  classification 
process.  These  are  body  waves,  travelling  through  the  earth,  and 
surface  waves,  travelling  along  the  earth's  surface.  Since  both  earth- 
quakes and  nuclear  explosions  generate  short-period  body  waves  and 
long-period  surface  waves,  the  conventional  criteria  used  to  distinguish 
between  earthquake  and  explosion  sources,  such  as  the  widely  used 

surface-to-body  magnitude  ratio  (M  : m ) criterion,  depend  on  freeing 

s b 

these  surface  and  body  waves  as  much  as  possible  from  interfering 

noise  and  other  waves  from  overlapping  events. 

The  detection  threshold  for  events  has  been  reduced  to  body 

magnitudes  (m  ) on  the  order  of  4.  0 by  the  advent  of  large  aperture 
b 

arrays  such  as  LASA,  NORSAR  and  ALPA.  Similarly,  the  development 
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of  very-low  period  high-gain  seismometers  promises  similar  reductions 
for  detecting  surface  waves  and  estimating  their  magnitude. 

One  of  the  major  problems  still  remaining,  is  estimating  the 
surface  wave  magnitude  (Mg)  of  an  event  when  overlapping  or  mixed 
events  are  present.  Since  mixed  events  are  often  a natural  occurrence 
for  multiple  earthquakes,  underground  testing  evasion  tactics  might 
include  timing  the  test  to  coincide  with  a local  earthquake  or  with 
the  coda  of  a large  earthquake.  According  to  the  Conference  of  the 
Committee  of  Disarmament,  on  the  average,  about  16  percent  of  the 
total  possible  single  station  observations  were  mixed  events,  in  which 
interfering  signals  from  an  overlapping  event  made  it  impossible  to 
extract  reliable  amplitude  or  spectral  information  from  the  wave  form. 

The  mixed  event  problem  can  be  more  clearly  defined  as 
separating  -the  Rayleigh  or  Love  primary  surface  waves  of  c. 
desired  event  from  a waveform  which  is  the  sum  of  the  desired 
event  and: 

1.  a Rayleigh  or  Love  (or  both)  primary  surface  wave 
of  an  interfering  event, 

2.  the  coda  of  a large  earthquake. 

Due  to  the  high  signal-to-noise  ratio  made  possible  by  the  VLP 

high  gain  seismometers,  OAS  proposed  to  investigate  signal  processing 
methods  which  were  capable  of  discrimination  of  mixed  events  in  a 
low  noise  environment.  Among  these  methods  were  included: 
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• it 

1. 

Polarization  Filter 

• • 

2. 

Azimuthal  Filter 

3. 

Homomorphic  Decompostion/Complex  Cepstrum 

« « 

4. 

Complex  Demodulation 

• • 

5. 

Matched  Filters 

6. 

Large  Array  Processing 

In  addition  to  these,  two  other  techniques  by  Alexander  concerning 
separation  of  Love  waves  from  Rayleigh  wave  interference  were 
found. 

As  an  extension  to  the  methods  of  Alexander  and  the  concepts 
of  the  polarization  and  azimuthal  filters,  OAS  has  developed  a filter 
termed  the  "ellipticity"  filter  which  has  the  capability  of  separating 
two  superposed  Rayleigh  waves  or  one  Rayleigh  and  one  Love  wave. 
Briefly,  the  ellipticity  filter  processes  the  three  signal  outputs  of  a 
three  component  seismometer  station  using  the  value  of  ellipticity  at 
the  station  to  produce  the  complex  amplitudes  and  azimuths  of  the 
two  superposed  surface  waves. 

The  record  is  divided  into  a number  of  intervals,  the  interval 
length  being  determined  by  the  time  variability  of  the  interfering 
wave.  For  each  interval,  the  waveform  is  transformed  into  the 
frequency  domain.  Each  frequency  bin  is  now  defined  by  three  complex 
numbers,  one  each  for  the  vertical  and  the  two  horizontal  channels. 
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Using  the  value  of  ellipticity  of  the  station  site,  the  superposition 
equations  for  the  case  of  two  Rayleigh  waves  or  one  Rayleigh  and 
one  Lovo  wave  may  be  written  for  each  frequency  bin.  These 
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equations  have  been  solved  in  closed  form  giving  the  two  complex 
source  amplitudes  and  their  respective  azimuths  in  terms  of  the 
three  complex  data  points  and  the  ellipticity. 

Having  solved  the  equations  for  all  frequency  bins,  the 
inverse  transform  of  those  complex  amplitudes  associated  with  one 
source  azimuth  value  yields  the  surface  waveform  for  that  source. 
Similarly,  the  inverse  transform  of  the  complex  amplitudes  associ- 
ated with  the  other  source  azimuth  yields  the  second  surface 
waveform. 

Similarly  derived  waveforms  for  ocher  time  intervals  in  the 
record  may  be  combined  smoothly  to  yield  the  overall  surface 
waveforms. 

The  main  limitations  in  this  method  are  those  imposed  by  the 
additions  of  ambient  noise  to  the  superposed  surface  waves.  The 
superposition  equations  have  been  solved  in  a deterministic  manner 
and  additive  noise  will  clearly  alter  the  accuracy  of  the  results. 

In  addition,  the  noise  sensitivity  of  the  solution  may  be  greater  for 
small  azimuthal  source  separations  than  for  large  separation.  These 
and  other  aspects  concerning  additive  noise  are  being  investigated 
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to  yield  limitations  of  this  approach.  This  investigation  includes 
both  theoretical  studies  and  empirical  studies  and  simulations. 

In  addition,  the  method  is  based  on  the  fact  that  the  ellipticity 
remains  constant  for  all  azimuth  around  the  station  site.  Nothing  in 
the  literature  has  been  found  to  prove  or  disapprove  this  assumption 
for  the  VLP  waves  in  question.  Thus,  we  are  initiatiig  an  effort  to 
determine  whether  or  not  ellipticity  exhibits  azimuthal  dependence  at 
a given  site  by  calculating  the  ellipticity  for  a number  of  earthquakes 
of  differing  azimuths  at  the  site. 

Section  2 of  this  report  reviews  the  methods  investigated  to  date, 
namely  the  polarization  and  azimuthal  filters  and  the  two  methods  by 
Alexander.  The  other  signal  processing  methods  listed  are  being 
looked  into.  Section  3 formulates  the  problem  mathematically  and 
gives  the  solution  method  comprising  the  ellipticity  filter  approach. 
Section  4 describes  some  preliminary  computer  simulation  testing  of 
the  filter  while  Section  5 outlines  areas  for  future  investigation  and 
testing  of  the  filter. 
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for  the  VLP  waves  in  question.  Thus,  we  are  initiating  an  effort  to 
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Review  of  Existing  Techniques 


The  ellipticity  filter  concept  is  an  extension  of  the  azimuthal 
and  polarization  filter  techniques  of  Choy  and  McCamy4as  well  as 
waveform  separation  algorithms  presented  by  Alexander.  These 
techniques,  in  addition  to  others,  are  reviewed  in  the  following 
section  with  emphasis  beirg  placed  on  their  utilization  in  the 
mixed  event  problem. 

2.  1 The  Polarization  Filter 

The  polarization  filter,  as  developed  by  Choy  and 
McCamy,  was  designed  to  separate  a 90°  polarized  Rayleigh  wave  from 
interfering  waves  by  retaining  information  in  those  time-frequency  bins 
which  possess  the  desired  attribute  of  a 90°  phase  shift  between  the 
vertical  and  ,rizontal  components  of  a three  component  seismometer. 
This  filter,  thus,  utilizes  those  advantageous  moments  in  which  the 
interference  is  small  compared  to  the  signal. 

For  a given  time  interval,  the  filter  computes  the  Fourier 
transform  of  the  vertical  and  one  of  the  horizontal  components  of  a 
three-component  seismometer.  For  each  frequency  component  in  the 
desired  passband,  the  polarization  filter  generates  a filter  function 
proportional  to  the  phase  difference  between  the  vertical  and  the 
selected  horizontal  component.  The  filter  function  used  by  Choy  and 
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McCamy  is 


F(uj)  = sinN[  - A (4j)l 

Z X 

where  uj  ) and  ^(*0  are  the  phases  of  the  vertical  and  selected 
horizontal  components  respectively  and  N is  a constant  which  deter- 
mines the  fall-off.  The  filtered  output  is  calculated  by  multiplying 
the  Fourier  transform  of  the  signal  by  F(oj).  Thus,  a component 
having  a phase  difference  near  90°,  as  for  a pure  Rayleigh  wave, 
will  be  passed,  whereas  non-Rayleigh  type  signal*  will  be  rejected. 

The  inverse  transform  of  the  filtered  components  is  taken  to 
yield  a time  waveform  for  the  interval.  This  procedure  is  repeated 
for  a number  of  overlapping  time  intervals  required  for  the  entire 
signal.  The  resultant  filtered  waveforms  are  added  using  a smoothing 
algorithm  to  yield  a continuous  time  function. 

The  main  disadvantage  of  this  approach  seems  to  be  that,  even 
a small  amount  of  interference  in  a.  given  time -frequency  bin  may 
sufficiently  corrupt  the  waveform  so  that  it  is  partially  rejected.  If 
two  Rayleigh  waves  are  superposed,  the  filter  will  either  attenuate 
or  pass  the  mixed  signal  depending  on  the  phase  difference  of  the 
mixed  signal.  In  either  case,  the  desired  event  has  not  been  separated. 
In  the  Rayleigh: Love  mixed  event,  the  Rayleigh  wave  will  pass  only  if 
it  significantly  dominates  the  Love  wave.  If  the  Love  wave  dominates, 
then  the  record  will  be  attenuated. 
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The  Azimuthal  Filter 


The  azimuthal  filter,  also  developed  by  Choy  and  McCamy, 
was  intended  to  separate  a Love  wave  from  interfering  waves  by 
retaining  information  in  those  time-frequency  bins  in  which  the 
azimuth,  calculated  from  the  two  horizontal  components,  is  near  to 
an  a priori  azimuth  value. 

The  same  procedure  is  followed  as  for  the  polarization  filter 
except  that  the  two  horizontal  components  are  processed.  A filter 
function  F(ca/)  is  generated  by  the  filter  and  has  the  form 

F(uli)  = sin  [ ob  - tan  X( ^ ) 

L X(uj) 

Thus,  if  the  azimuth  of  the  signal  component  is  near  that  expected, 
the  component  is  passed.  If  not,  the  component  is  attenuated. 

The  azimuthal  filter  is  intended  to  separate  Love  waves  from 
Rayleigh  by  passing  the  Love  wave  when  it  dominates  the  signal. 
However,  even  the  presence  of  a smaller  interfering  signal  can 
perturbate  the  azimuth,  thus  attenuating  the  record.  For  instance, 
the  particle  motion  of  a Rayleigh  wave  from  the  same  azimuth  as  the 
Love  wave  is  orthogonal  to  the  particle  motion  of  the  Love  wave.  Thus, 
if  the  Rayleigh  wave  has  an  amplitude  of  -6dB  with  respect  to  the  Love 
wave,  the  apparent  azimuth  of  the  combined  wave  is  27°  from  the 
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actual  azimuth  value  for  the  Love  wave.  A coefficient  of  N = 8 in  the 


filter  function,  F(v^),  causes  a 5dB  attenuation  of  this  component. 


2.  3 Frequency  Dependent  Rotation  of  Axis 


A method  suggested  by  Alexander  would  separate  Love  waves 


from  interfering  Rayleigh  waves  by  a frequency  dependent  rotation 


of  the  horizontal  axis.  This  method  requires  the  previously  deter- 


mined azimuth  of  arrival  of  the  Rayleigh  interference  on  a frequency 


dependent  basis.  The  azimuthal  dependence  on  frequency  is  based  on 


the  lateral  refraction  of  Rayleigh  waves  off  ocean-continen*-  boundaries 


being  frequency  dependent.  Using  this  azimuth  information,  a node  is 
placed  in  the  direction  of  the  azimuth  of  arrival  on  a per  frequency  basis 
by  a suitable  rotation  of  horizontal  coordinates.  This  node  eliminates  the 
Rayleigh  components  leaving  only  the  desired  Love  components  (assuming 
that  the  azimuth  of  arrival  of  the  Love  wave  is  not  90  from  the  Rayleigh). 


This  method  would  clearly  be  applicable  to  the  mixed  event  in 


which  two  Rayleigh  waves  are  superposed,  assuming  their  azimuths 


are  not  identical.  The  only  disadvantage  lies  in  determining  the 


azimuth  of  arrival  of  the  event  to  be  eliminated.  Alexander  suggests 


that  an  array  using  F-k  processing  would  be  a likely  method. 
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2.4  Frequency  Dependent  Subtraction 

of  Rayleigh  Interference  using  the  Ellipticity  Constraint 

A second  method  proposed  by  Alexander  for  separating  Love 
waves  from  Rayleigh  interference  is  to  subtract  out  the  Rayleigh 
components  in  the  horizontal  signals  using  the  azimuth  of  arrival  of 
the  Rayleigh  signal  and  the  ellipticity  constraint  between  the  vertical 
and  horizontal  components  of  a Rayleigh  wave.  Sin  = e the  vertical 
component  contains  only  the  Rayleigh  interference,  the  azimuth  of  arrival 
and  the  ellipticity  constraint  define  the  amount  of  Rayleigh  interference 
contained  in  each  horizontal  component.  When  the  interfering  components 
are  subtracted,  the  resulting  signal  should  be  the  desired  Love  com- 
ponents. This  procedure  would  also  be  done  on  a per  frequency 
basis  since  the  azimuth  of  arrival  of  the  Rayleigh  interference 
may  be  frequency  dependent. 

As  in  the  previous  method,  the  azimuth  of  the  interfering 
Rayleigh  wave  must  be  accurately  determined  to  yield  successful 
results.  In  addition,  the  ellipticity  constraint  at  the  station  site 

must  be  learned  in  order  to  subtract  out  the  correct  amount  of 
interference. 

Utilizing  this  method  for  the  Rayleigh:Rayleigh  mixed  event 
problem  implies  determination  of  the  complex  amplitude  and  the 
azimuth  of  the  interfering  Rayleigh  wave  since  the  vertical  signal  now 
contains  components  of  both  signals. 
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3.  0 The  Seismic  Ellipticity  Filter 


In  an  effort  to  utilize  the  information  represented  by  a three- 
component  seismometer  station,  a decomposition  algorithm  termed 
the  "Seismic  Ellipticity  Filter"  has  been  developed.  As  input,  this 
filter  accepts  the  superposition  of  two  independent  surface  waves 
(either  Rayleigh  or  one  Rayleigh  and  one  Love)  from  different 
azimuths.  Using  the  ellipticity  constraint  defining  the  radial-to 
vertical  amplitude  ratio  of  a Rayleigh  wave,  the  ellipticity  filter 
operates  on  this  input  and,  in  the  noiseless  case,  perfectly  separates 
the  superposed  signals  into  the  two  original  waveforms.  In  addition, 
the  azimuths  of  the  two  wavei  are  determined. 

In  general,  the  filter  operation  is  valid  for  high  signal-to- 
noise  ratios  and  initially  solutions  are  found  only  for  the  noise-free 
case. 

This  filter  is  an  outgrowth  of  the  azimuthal  and  polarization 
filter  concepts  of  Choy  and  McCamy  and  an  extension  of  the  method 
of  Alexander,  discussed  in  Section  2.4,  to  utilize  the  ellipticity 
constraint  in  separating  Love  waves  from  Rayleigh  interference.  In 
our  case,  the  ellipticity  constraint  is  used  to  formulate  the  superposition 
equations  defining  the  three  seismometer  signals.  This  formulation 
allows  us  to  solve  both  the  Rayleigh  wave  with  superposed  Rayleigh 
wave  interference  problem  as  well  as  the  Love  wave  with  superposed 
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I 

Rayleigh  wave  interference  problem. 

The  desired  wave  can  either  be  the  Love  or  Rayleigh,  while 
the  interfering  wave  is  generally  considered  Rayleigh.  The  most 

« • 

frequent  problem  is  that  of  an  event  being  masked  by  earthquake  coda 
which  is  considered  to  be  Rayleigh  in  character.  In  the  following 
development,  the  two  waves  will  be  designated  simply  by  Source  1 
and  Source  2,  since  the  filter  solves  for  both  waves  and  does  not 
care  which  is  signal  and  which  is  interference. 

3.  1 Problem  Statement 

Given  the  outputs  of  three  seismometer  components  (E/W, 

N/S,  and  vertical),  devise  a method  which  can  extract  a 
desired  primary  surface  wave,  either  Rayleigh  or  Love, 
from  signals  which  are  the  superposition  of  the  desired 
wave,  with  an  interfering  Rayleigh  surface  wave. 

Since  the  filter  operates  in  the  frequency  domain,  the  three 

seismometer  components,  E/W,  N/S,  and  vertical,  are  represented 

by  their  Fourier  Transforms  as  X,  Y,  and  Z respectively,  and  X 

denotes  the  j**1  frequency  component  of  X. 

The  following  assumptions  concerning  signal  and  filter 

parameters  are  made: 

1.  Only  two  unknown  signals  are  present  in  any  specific 

time  interval  of  observables  being  processed.  These  signals 
must  either  be  both  Rayleigh  or  one  Rayleigh  and  one  I. eve. 

The  superposition  equations  which  define  the  observable 

quantities  are  limited  by  the  number  of  independent  seismometer 

components,  in  this  case,  three.  Each  of  the  observables,  X.,  Y,, 
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I 

I 

I 

I 


I 

I 

I 

I 
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and  Zj  is  a complex  quantity  and,  therefore,  represents  two  real 
superpr  sition  equations.  Altogether,  there  are  six  real  superposition 
equations  which  define  the  three  complex  observables. 

Each  unknown  signal  represents  three  real  unknowns:  an 

amplitude,  phase  and  azimuth.  Thus,  two  unknown  signals  represent 
six  real  unknowns,  the  same  number  as  the  number  of  real  super- 
position equations.  Hence,  more  than  two  unknown  signals  could  not 
be  solved  by  using  only  three  independent  seismometer  components. 
Similarly,  two  Love  waves  cannot  be  solved  as  the  problem  contains 
six  unknowns  but  only  four  real  equations,  since  the  vertical  component 
contains  no  information. 

Due  to  the  dispersive  nature  of  seismic  surface  wave  signals, 
for  a time  interval  small  compared  to  the  primary  surface  wave 
duration,  most  of  the  energy  of  the  surface  wave  in  that  interval  will 
be  concentrated  in  a few  spectral  components.  Thus,  a third  unknown 
signal  may  possibly  be  resolved  if  it  is  significantly  shifted  in  time 
such  that  its  contribution  to  the  major  frequency  components  for  the 
two  unknown  signals  is  negligible. 

2.  The  observables  arc  noise-free. 

The  superposition  equations  defining  the  ellipticity  filter  do  not 
contain  explicit  noise  terms  but  are  a set  of  six  deterministic  equations 
in  six  unknowns.  Clearly  the  closed  form  solutions  will  be  degraded 

in  the  presence  of  noise.  It  is  felt,  though,  that  a high  signal -to- 
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noiee  ratio  environment  will  still  yield  satisfactory  resi  Its.  The 

threshold  for  such  operations  has  yet  to  be  determined. 

3.  The  apparent  ellipticity  at  the  seismometer  station 

site  is  known  and  does  not  vary  with  azimuth  of  arrival. 

The  apparent  ellipticity  is  a complex  quantity  which  equals  the 
radial  component  divided  by  the  vertical  component  of  a pure  Rayleigh 
wave,  where  these  quantities  are  measured  at  the  output  of  the  seis- 
mometers. Thus,  any  seismometer  phase  shift  or  amplitude  scale 
factor  is  included  in  this  factor. 


Seismic  Filter  Configuration 


The  seismic  filter  conifguration  is  shown  in  Figure  3.  1. 
This  pictorial  representation  shows  the  method  followed  for  the  k* 
time  interval. 


The  three  seismometer  component  waveforr.  s,  x(t),  y(t),  and 
z(t),  are  sampled  in  time  to  yield  x^(t),  y^(t),  and  z.(t).  These  wave- 
forms are  divided  into  either  adjacent  or  overlapping  time  intervals, 
the  k**1  interval  denoted  by  a superscript  k.  The  duration  of  this 

time  interval  depends  mainly  on  the  time  variability  expected  in  the 

3,4 

azimuth  of  arrival.  For  coda-type  interference,  the  consensus 

appears  to  be  intervals  of  approximately  2-3  minutes. 

For  the  time  intervals  selected,  the  Fast  Fourier  Transforms 

k k k k 

of  x , y , and  z are  taken  to  yield  X , Y , and  Z.  , where  X. 
i i i J j J J 
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denotes  the  frequency  component  in  the  k time  interval  of  x(t). 


For  each  frequency  component  in  each  time  interval,  solve  the 


superposition  equations  (see  Sections  3.3  and  3.4)  for  the  two  unknown 


complex  amplitudes  A^  and  and  their  respective  azimuths  0^  and 


In  the  k**1  interval  associated  with  Source  1,  are  those  complex 


k 

amplitudes  having  azimuth  0^  . Associate  with  Source  2 all  those 


k 

complex  amplitudes  having  azimuth  0^  . Take  the  inverse  Fourier 


Transform  of  those  complex  amplitudes  associated  with  Source  1, 


yielding  the  signal  a^  (t).  Similarly,  take  the  inverse  Transform 
of  all  those  complex  amplitudes  associated  with  Source  2,  yielding 


*>• 


Over  all  k of  interest,  smooth  those  resultant  time  waveforms 


associated  with  Source  1 in  which  0^  is  the  same.  The  resultant 


waveform  is  the  total  waveform  for  Source  1.  Repeat  for  Source  2. 


3.3  The  Seismic  Ellipticity  Filter  for  Rayleigh:  Rayleigh  Mixed  Events 


Based  on  the  source  geometry  given  in  Figure  3.2,  the  super- 

k k k 

position  equations  defining  the  observables  X.  , Yj  , and  Zj  for  two 


Rayleigh  sources  are  given  by 
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k k k 

X.  = -f(  A.jj  sin  0j  + A2j  sin  02  ) 


Y. 

J 


-f(  Ajj  cos  0 j + A^j  c°s02  ) 


k k 

A . . + A 

lj  2j 


where  f is  the  apparent  ellipti.city  at  the  station  site. 


These  three  complex  equations  actually  represent  six  equations  by 

separating  real  and  imagina-v  parts.  There  are  six  unknowns: 
k k 

Ajj  and  A^  are  complex  and,  thus,  represent  two  unknowns  each 

(an  amplitude  and  a phase);  the  two  real  azimuths  0 and  0 . Thus. 

1 2 

we  have  six  equations  with  six  unknowns.  These  equations  have  been 

solved  (See  Appendix  A)  in  closed  form.  For  simplicity,  we  drop 

the  subscripts  j and  k.  The  complex  amplitudes  A ^ and  A^  are  given  by 


A2  = |a2I  e1^2 

Thus,  the  six  unknowns  , A2  ' ^1  » ' 61  ' and  82  are  given  bY 

Im(Z)  cosf>2  - Re(Z)  sind2 

®in(  - tz) 


I 


North 


Source  1 


Amplitude  or 


East 


Source  2 


Amplitude  A^  or 


where  6j,  02  are  the  real  azimuths  of  the  two  sources, 


Aj,  A2  are  the  complex  vertical-component  amplitudes  for 
LR  (Rayleigh)  sources 


Bj,  B2  are  the  total  complex  horizontal  amplitudes  for  LQ 

(Love)  sources  - positive  sense  to  the  right  of  phase 
velocity. 


Fig.  3.2:  Source  Geometry 


J 


where 


set  of 
tities  I 
source 


a2!  ■ 


Im(Z)cos$i  - Re(Z)sin($i 
sin(  ) 


*1 


-b  + (b  - 4ac) 


1/2 


2a 


-b  - (b*  - 4ac) 


1/2 


2a 


arc  tan 


Im(X)  - P,,  Re(X) 

£« 

Im(Y)  - P Re(Y) 
2 


0 = arc  tan 

2 


Im(X)  - Pj  Re(X) 
Im(Y)  - P1  Rc(Y) 


arc  tan  t)  , 
rl 


P = arc  tan  i 
2 2 


r- Y 


X 


= Re  (y)  + Re  (y)  - Re  (Z) 


b = 2[  Re(Z)  Im(Z ) - Re(  y ) Im(  y-  ) - Re(  y- )Im(y ) ] 


2 Y 2 x 2 

Im  ( — ) + Im  ( -£■ ) - Im  (Z) 


Thus,  for  each  frequency  band  in  each  time  interval,  the  above 

squations  is  solved  for  the  two  sources  and  yields  the  six  quar.- 

A I , 0 , and  0 for  source  1 and  i A | , i,  and  0 for 
111  2 Z b 

2.  The  solutions  for  all  frequency  bins  in  a time  interval  are 
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combined  by  associating  with  source  1 all  those  solutions  having 

azimuth  6 . Similarly,  those  solutions  having  azimuth  0 are 
1 2 

associated  with  source  2.  The  time  waveform  in  the  time 
interval  for  source  1 is  then  calculated  by  summing  those  sinusoids 
associated  with  source  1.  Thus, 

Si  (t)  = ElAj  sin  v ui.  t + <j)  ) 

1 ? 3 lj 


Similarly, 


(t) 


■Ei 


2j 


sin(  uj.  t + 
J 


j k 

*Zj  > 


This  procedure  is  applied  over  all  time  intervals  defining  the  record. 

The  overall  signals  s (t),  [or  s (t)  ],  are  then  formed  by  summing 

1 £ 

k 

in  some  smooth  manner  the  interval  waveforms  s (t)  over  k. 


3. 4 The  Seismic  Ellipticity  Filter  fc>*  Raleigh:  Love  Mixed  Events 


For  the  Love:Rayleigh  mixed  event,  the  superpositior  equations 

k k k 

defining  the  observables  X.  , Y.  , and  Z;  are  given  by 

J 1 J 


■f  A sin0 

1 


B cos  0 
2j  2 


Y.  = -f  A k cos  0 + B sin  0 

J lj  1 2j  2 


Zk  . Ak 

j lj 
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where  f is  the  apparent  ellipticity  and  represents  the  complex 
amplitude  of  the  Love  wave. 

As  in  Appendix  B,  these  three  complex  equations  represent  six  real 
equations  in  the  six  unknowns  Re  A^,  Im  A^,  Re  B^,  Im  B^,  0^  and  0^ 


0_  = arc  tan 

2 


Im(Y)  - P»Re(Y) 
Im(X)  - P.Re(X)  J 


Im  B 


Im  B^  = 


0 + arc  cos 

2 


sin0 


sin0. 


[ 

M 

['"'r 


Re(  y*)  sin02  + Re(-j)  cos  02 


Re(Z) 

) + Re(Z)  cos  0^ 


) + Im(Z)  cos  0 


1 


where 


ReA  j = 

tmAj  = 


Im(Z) 

Re(Z) 

Re(Z) 

Im(Z) 


Combining  the  solution  for  all  frequency  bins  and  time  segments 
to  yield  the  source  waveforms  is  performed  as  outlined  in  Section  3.3 
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Some  preliminary  analysis  of  the  ellipticity  filter's  operation 
was  performed  using  simulated  input  data.  A time  interval  of  180 
seconds  was  chosen.  Using  that  interval,  the  3rd  through  7th  harmonics 
(29.9  sec,  30.0  sec,  36.0  sec,  45.0  sec,  60.0  sec  periods)  were  used 
to  construct  the  Airy  phase  of  a Rayleigh  waveform.  In  addition,  an 
interfering  coda  waveform  was  similarly  constructed.  The  computer 
program  mixed  these  two  waveforms  assuming  their  propagation  paths 
were  orthogonal.  The  Signal,  Coda,  and  mixed  Signal  + Coda  are 
shown  in  Figures  4.  la,  b,  and  c.  These  are  pictures  of  the  vertical 
component.  The  horizontal  components  are  not  shown.  The  computer 
program  perfectly  separated  the  mixed  Signal  + Coda  in  the  Airy  signals 
labeled  Output  1 and  Output  2,  as  in  Figure  4.  Id  and  e.  Note  that 
Output  1 is.  identical  to  the  input  signal  while  Output  2 is  identical  to 


the  coda. 

The  actual  computer  output  is  shown  in  Figure  4 . 2.  This  shows 
seven  separate  computer  runs  for  seven  different  signal-to-coda  ratio 
from  -9  dB  to  +9  dB  in  3 dB  steps.  Each  run  lists  the  azimuth  and 
the  sine  and  cosine  components  of  the  complex  amplitude  for  each  of 


22 


roftFO 


ccrA 

.3536 

Kri«lon 

bO.O 

bE.C 

45.0 

SEC 

noise 

. Olbo 

tc  Ii' UH 

‘>0.0 

-.0  0 

-.0 

9o.no 

Mi.  cOH 

• 00 

.01 

.07 

.35 

CCS  COP 

-05 

1.01 

-.bo 

.00 

COOA 

. be  HO 

rE* IOD 

bO  . 0 

SEC 

45.0 

SF.C 

NOISE 

.0156 

A 2 I MM  H 

90 . n 

-.0  J 

-.0 

90. OC 

SI).  uC.vp 

•01 

.01 

• 07 

.50 

coj  ocvr 

-.40 

1.01 

-.bC 

.00 

CODA 

.71)71 

HErtI JD 

60. 0 

SEC 

4b. 0 

SEC 

NOISE 

.Olbo 

A2IPI  IP 

90. n 

-.00 

.0 

90.00 

SIN  CONP 

• 01 

.01 

. d7 

.71 

CCb  LC1* P 

-.70 

1.01 

-.bU 

.01 

COl  a' 

1.0000 

PCH I uO 

bO.O 

SEC 

43.0 

SEC 

NOISE 

• 0156 

A2 I VL  1 H 

. 0 

-.00 

-.0 

90.00 

SIu  CCNP 

•01 

.01 

.07 

1.00 

COb  CO*P 

-.99 

1.01 

-•bO 

.01 

CODA 

1.4142 

HTHloO 

oO.  0 

SEC 

45.0 

SEC 

NCISE 

. 0 1 bb 

A 2 li»'  0 I H 

90.0 

- • 0 J 

.0 

90.00 

SIN  LCMP 

• 01 

.01 

.e7 

1.41 

COS  CO«P 

-1.40 

1.01 

-.50 

• Ol 

COCA 

2.000U 

KE  H I JD 

t-0.0 

SEc 

45.  C 

SEC 

NOISE 

• Olbo 

A 2 I * U 1 H 

90. C 

-.00 

-.0 

90.00 

sir.  v.ovp 

• 02 

.01 

.07 

2.00 

CCb  OOKP 

-1 .9P 

1.01 

-.50 

.02 

CODA 

2.0204 

f'EhlOC 

60.0 

SEC 

45 . <1 

SEC 

NOISE 

.Olbo 

A^lNUlH 

9U.0 

-.0  3 

-.0 

90.00 

SIN  v.Ci*P 

• 03 

.01 

.07 

2.P3 

COb  cCi'-P 

-200 

1.01 

-.50 

.03 

.FI'. 


Kip.  4.2: 


i3 


45.0 

SEC 

36.0 

SEC 

3(1.0 

-.0 

90.  nn 

98.0 

.00 

90.0 

• d7 

.35 

• 00 

.50 

.35 

•.5u 

.00 

-•3b 

-.07 

.00 

45.0 

SEC 

36.0 

SFC 

30,0 

-.u 

90.00 

90.0 

.00 

9U.0 

. a7 

.50 

.00 

.50 

.50 

•.50 

.00 

-.50 

-.67 

.00 

45.0 

SEC 

36.0 

SEC 

30.0 

.0 

90.00 

90.0 

.00 

90.0 

.37 

.71 

.01 

.50 

.71 

-.50 

.01 

-.71 

-.67 

.01 

46.0 

SEC 

v 36 . 0 

SFC 

30.0 

-.0 

90.00 

90.0 

.00 

90.0 

.ti7 

1.00 

.01 

.50 

1.00 

-.50 

.01 

-1.00 

-.87 

.01 

45.0 

SEC 

36.0 

SEC 

30.0 

.0 

90.00 

90.0 

.00 

90.0 

.e7 

1 .4 1 

.01 

.50 

1.41 

.50 

• 01 

-1.41 

-.87 

.01 

45. C 

SEC 

36. 0 

SEC 

30.0 

-.0 

90.00 

90.0 

.00 

90.0 

.87 

2.00 

.02 

.50 

2.00 

.50 

.02 

-2.00 

-.87 

.02 

45.0 

SEC 

36.0 

SEC 

30.0 

-.0 

90.00 

90.  0 

.00 

90.0 

.87 

2.83 

.03 

.50 

2.63 

.50 

.03 

-2.d3 

-.87 

.03 

SEC 

25. 7 

SFC 

.09 

-►0 

90.00 

.01 

-.67 

.00 

1.00 

-.50 

• 

-.35 

SFC 

25.7 

SFC 

.00 

-.0 

90.00 

.01 

-.07 

.00 

l.on 

-.50 

-.50 

SEC 

25.7 

SEC 

.00 

-.0 

90.00 

.01 

-.87 

.01 

1.00 

-.50 

-.71 

SEC 

25.7 

SEC 

.00 

-.0 

90.00 

.d 

-.87 

.01 

1.00 

-.50 

-1.00 

SEC 

25.7 

SEC 

.00 

-.0 

90.00 

.01 

-.87 

.01 

1.00 

-.50 

-1.41 

SEC 

25.7 

SFC 

.00 

-.0 

90.00 

.01 

-.e7 

.0? 

1.00 

-.50 

-2.00 

SEC 

25.7 

SFC 

.00 

-.0 

90.00 

.01 

-.87 

.0* 

1.00 

-.50 

-2.8? 

Reproduced  from 
best  available  copy. 


Fig.  4.2:  Computed  Filter  Output 
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Future  Investigation  of  the  Ellipticity  Filter 


Performance  evaluation  of  the  ellipticity  filter  involves  testing 

• i 

in  a number  of  areas.  One  of  the  most  important  is  the  inclusion  of 

I 

additive  noise  in  the  mixed  event  waveforms  to  determine  the  noise 

level  which  can  be  handled  before  the  filter  algorithm  breaks  down. 

I 

Secondly,  the  filter  may  be  more  susceptible  to  additive  noise  when 
the  azimuth  of  the  mixed  events  are  relatively  near.  Thus,  azimuth 
of  the  two  sources  will  be  a parameter  in  the  noise  investigation.  A 
third  factor  may  be  the  relative  strength  of  the  two  signals  in  the 
presence  of  additive  noise. 

Successful  filter  operation  has  been  based  on  exact  knowledge 
of  the  ellipticity.  When  processing  real  data,  the  ellipticity  value  used 
may  differ  from  the  true  value.  Thus,  we  will  investigate  the  effect  of 

variations  in  the  value  of  ellipticity  used  on  filter  operation.  In  addition, 

. 

the  value  of  ellipticity  at  a station  site  has  been  assumed  independent 
of  the  azimuth  of  arrival.  A test  will  be  performed  to  determine  the 
actual  value  of  ellipticity  at  a station  site  and  its  sensitivity,  if  any,  to 
the  azimuth  of  arrival. 

The  above  investigations  will  be  performed  using  both  simulated 
and  real  VLP  high- gain  data  from  the  Lamont- Doherty  Geological  Observ- 
atory of  Columbia  University.  Presently,  we  are  compounding  a data 
bank  of  explosion  events,  earthquakes,  earthquake  coda  and  ambient  noise 
for  testing  purposes. 
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APPENDIX  A : Derivation  of  the  Ellipticity  Filter 

for  Two  Mixed  Rayleigh  Events 


Assume  we  have  a three  component  seismometer  with  the  output 
signals  x{t),  y(t),  and  z(t>  being  the  E/W,  N/S  and  vertical  components 
respectively.  Assume  that  a Fourier  Transform  of  these  signals  has 
been  taken  for  some  desired  time  window.  The  signals  X,  Y,  and  Z 
represent  one  of  the  Fourier  components  for  the  tr- informed  signals 
and,  thus,  are  complex  functions  of  frequency. 

For  the  case  of  two  mixed  Rayleigh  sources,  the  superposition 
equations  for  X,  Y,  and  Z may  be  written  as 


where 


X = 

-f  ( 

Ai 

sin  0j 

+ A,  sin  0 ) 

2 2 

(1) 

Y = 

-f  ( 

Ai 

COS  0 j 

+ A 2 cos  02  ) 

(2) 

Z = 

A1 

+ 

A2 

(3) 

A1  and  are  the  complex  amplitudes  of  the  two 
Rayleigh  events, 

01  and  02  are  the  azimuths  of  Aj  and  A2,  respectively, 
f is  the  apparent  ellipticity. 


The  source  geometry  is  given  in  Figure  3.  1,  where  f is  the  apparent 
ellipticity  at  the  receiving  site  and  is  assumed  know  j. 

From  (2) 

Y 

Re  ( - - ) = Re  ( Ax  ) cos  0J  + Re  ( A2  ) cos  02  ( 

Y 

Im  ( - — ) = Im  ( A ) cos  0 + Im  ( A ) cos  0 ( 

1 1 2 2' 
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a 


- 


Im  ( A i ) 


yields 


Multiplying  (4)  by  = 

Re  (Aj  ) 

Y 

pl  • Re  ( - j ) = Im  (AJ  )cos  + P ^ Re  (A2 ) cos  (6) 

Subtracting  (6)  from  (5) 

(-  J-  ) - Pj'Re  ( - J-  ) = [ Im(A2)  - P-Re  (A2)  ] cos  02  (7) 

Substituting  into  (7)  the  value  of  A from  (3), 

c* 

Im(  - y ) - P^  Re(  - y ) = [ Im ( Z ) - P^RetZ)  ] cos  ©2  (8) 

By  a similar  sequence,  equation  ( 1)  yields: 

Im(  - y ) - Pj*  Re  ( ~ ) = [ Im(Z)  - P^  Re  (Z)  ] sin  ©2  (9) 

Squaring  (8)  and  (9)  and  adding  gives: 

Im2(|-)  - ZPjImt.Y)  Re  ( Y_  ) + p2  Re2  (' y ) 

2 X X v 2 2 v 

+ Im  (—  ) - 2PJ  Im  (—  ) Re  (-p)  + Pj  Re  ( ~ ) 

= [P2Re2(Z)  - 2P  Re(Z)  Im(Z)  + Im2(Z)  ] (cos  0 + sin  0 ) 

11  2 2 

(10) 

2 2 2 2 
R®  ( ) = [Re(  )]  and  Im  ( ) = [lm(  )] 

(10)  yields  the  quadratic  form: 

2 

aPl  + bPl  + C = 0 (11) 


where 


Rearranging 
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where 


2 Y 2 X 2 

= Re  ( y)  + Re  ( — ) - Re  ( Z ) 


b = 2[  Re  ( Z ) Im  ( Z ) - Re  (I  - Re  ( — )Im(£)l 

f f f f J 

c = Im2(y)  + Im2(  j)  - Im2(Z) 


By  an  entirely  similar  derivation,  P , defined  by 

2 


P = 

2 Re  (Z2  ) 

satisfies  the  same  quadratic  equation. 

Thus,  and  P^  are  the  two  roots  of  the  same  quadratic  equation. 
We  arbitrarily  assign  the  two  roots  as 


2 1/2 

-b  + ( b - 4ac  ) 


2 1/2 

-b  - ( b - 4ac  ) 


Dividing  (9)  by  (8): 


= arc  tan 


Im  (jr  ) 

Im  ( y"  ) 


Pi  Re(ir) 

Pj  Im(j) 


Similarily, 


= arc  tan 


X 

Im  (J) 
Im  (j-  ) 


p2  Re  (T) 


p2  Re  (TJ 


If  we  now  let 


, = 'Axi  •**  • -s 


then,  clearly. 


/j  = arc  tan  Pj  and  ^ = arc  tan  ^2 


Equation  (3)  becomes 


. U,l  ^ ♦ IaJ  A 


Multiplying  (16)  by  e and  taking  the  imaginary  part  of  the  result: 


-Re(Z)sin^2  + Im(Z)cos$2  = I sin(  ^ ^ ) 


and  | A \ 


Im(Z)cos^2  - Re(Z)sin(i 
sin  ( - 1>2  ) 


, Im(Z)  - P7Re(Z) 

-2  or  1 (17) 


sin^  - P2cos«(1 


Similarily 


Im(Z)cos$.  - Re(Z)sin^j  Im(Z)  - P Re(Z) 

or 


sin  ( 02  - ^ ) 


sin^  - P^cos£f2 


7 


APPENDIX  B: 


Derivation  of  the  Ellipticity  Filter 
for  Mixed  Rayleigh  and  Love  Events 


As  in  the  previous  derivation,  the  superposition  equation 
defining  the  observables  may  be  written  as 


X = 

-f  A j sin  0 

- B2  cos  02 

(1) 

Y - 

-f  A^  cos0 

+ Vin62 

(2) 

Z = 

A 

(3) 

1 

where  A is  the  complex  amplitude  of  the  Rayleigh  event, 

is  the  complex  amplitude  of  the  Love  event, 

0 , 0 are  the  azimuths  of  A and  B respectively, 
12  12 

f is  the  apparent  ellipticity. 

From  (1)  and  (2) 


Re(Y) 

= Re  ( -fZ  )cos  0 ^ + 

Re  (B  ) sin  0 
2 2 

(4) 

Im(Y) 

= Im  ( -fZ  )cos  0 + 

Im  (B^)  8*n  ®2 

(5) 

Let  P = Im(-fZ)/Re(-fZ).  Multiply  (4)  by  P: 

P-Re(Y)  = Im  ( -fZ)cos0,  + P-Re(B  J sin  0 (6) 

1 2 2 

Subtracting  (6)  from  (5): 

Im(Y)  - P*  Re(  Y)  = [ Im(B  ) - P Re(B  )]  sin  0 (7) 

m b fa 
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F rom  ( 3 ) 


Re  (X)  = 

Re  ( -fZ)  sin  0^ 

Re  (B  )cos  0 
2 2 

(8) 

Im  (X)  = 

Im  ( -fZ)sin  0 ^ 

Im  (B  )cos  0 
2 2 

(9) 

Multiplying  (8)  by  P and  subtracting  from  (9): 


Im(X)  - P*Re(X)  = [P  Re  (B  ) - Im(B  )1  cos  0 (10) 

2 2 2 

From  the  ratio  of  (7)  to  (10): 

e2  -Im(Y)  - PRe<Y>  (11) 

Im  (X)  - Re(X) 

Rearranging  (4)  and  (8)  and  taking  their  ratio: 


or 


or 


Re(Y)  + Re(fZ)cos0i 

2 = 1 (12) 

Re(X)  + Re(fZ)sin0^ 


Re(fZ ) [ cos 0 j + tan  0^  sinOj  ] = -Re(X)  tan02  - Re(Y)  (13) 

cos(  0 -0  ) = - Re(X)tar.02  4-  Re(Y) 

1 2 (14) 

sec  0 • Re  (fZ) 

2 


and  finally 


0 

1 


0 

2 


+ arc  cos 


Re(X)sin02  + Re(Y)cos02 
Re(fZ) 


(15) 
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From  (5) 


and  (1)  yields 


Im(Y)  + Im(fZ)cos0. 
Im  (B  ) l 

sin  0_ 


(16) 


Re  (B2)  = 

Re(Y) 

as 

Re  {Aj)  = 

Re  (Z) 

Im  (A  ) = 

« 

Im  (Z) 

sin  0. 


(17) 
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